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ABSTRACT 

To elucidate the effects of the UV background radiation on the collapse of pregalac- 
tic clouds, we implement a radiation-hydrodynamical calculation, combining one- 
dimensional spherical hydrodynamics with an accurate treatment of the radiative 
transfer of ionizing photons. Both absorption and scattering of UV photons are ex- 
plicitly taken into account. It turns out that a gas cloud contracting within the dark 
matter potential does not settle into hydrostatic equilibrium, but undergoes run-away 
collapse even under the presence of the external UV field. The cloud center is shown 
to become self-shielded against ionizing photons by radiative transfer effects before 
shrinking to the rotation barrier. Based on our simulation results, we further discuss 
the possibility of H2 cooling and subsequent star formation in a run-away collapsing 
core. The present results are closely relevant to the survival of subgalactic Population 
HI objects as well as to metal injection into intergalactic space. 
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1 INTRODUCTION 

It is widely recognized that the Ultraviolet (UV) background 
radiation, inferred from the proximity effect of Lya absorp- 
tion lines in QSO spectra (e.g., Bajtlik, Duncan & Ostriker 
1988; Bechtold 1994; Giallongo et al. 1996), is likely to ex- 
ert a significant influence upon the collapse of pregalactic 
clouds and consequently the formation of galaxies. Several 
authors have argued that the formation of subgalactic ob- 
jects is suppressed via photoionization and photoheating 
caused by the UV background (Umemura & Ikeuchi 1984; 
Ikeuchi 1986; Rees 1986; Bond, Szalay & Silk 1988; Efs- 
tathiou 1992; Babul & Rees 1992; Zhang, Anninos & Nor- 
man 1995; Thoul & Weinberg 1996). Final states of the pho- 
toionized clouds, however, are still unclear, because most of 
these studies assume optically-thin media, which precludes 
us from the correct assessment of self-shielding against the 
external UV fields. Self-shielding not only plays an impor- 
tant role in the thermal and dynamical evolution, but is also 
essential for the formation of hydrogen molecules (H2) which 
control the star formation in metal poor environments such 
as primordial galaxies. Several attempts have recently been 
made to take account of the radiative transfer of ionizing 
photons, adopting for instance a pure absorption approxi- 
mation (Kepner, Babul & Spergel 1997; Kitayama & Ikeuchi 
2000, hereafter KI), the photon conservation method (Abel, 



Norman & Madau 1999), and the full radiative transfer 
treatment (Tajiri & Umemura 1998, hereafter TU; Barkana 
& Loeb 1999; Susa & Umemura 2000). As for the dynami- 
cal states, Kepner, Babul & Spergel (1997) and Barkana & 
Loeb (1999) have considered hydrostatic equilibria of spheri- 
cal clouds within virialized dark halos. KI have explored the 
hydrodynamical evolution of a spherical system composed of 
dark matter and baryons. The full radiative transfer treat- 
ment, however, has not hitherto been incorporated with the 
hydrodynamics of spherical pregalactic clouds. 

In this paper, we attempt to implement an accurate 
radiation-hydrodynamical (RHD) calculation on the evolu- 
tion of spherical clouds exposed to the UV background, solv- 
ing simultaneously the radiative transfer of photons and the 
gas hydrodynamics. Our goals are to predict at a physically 
reliable level final states of photoionized clouds and also to 
assess accurately self-shielding against the UV background 
that is essential for H2 cooling and the subsequent star for- 
mation. Throughout the present paper, we assume the den- 
sity parameter Qq = 1, the Hubble constant h = Ho /(100 



km s Mpc 
n b = 0.1. 



0.5, and the baryon density parameter 
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2 MODEL 

A pregalactic cloud is supposed to be a mixture of baryonic 
gas and dark matter with the mass ratio of £lh '■ fio — = 
1:9. The numerical scheme for the spherical Lagrangian 
dynamics of two-component matter follows the method de- 
scribed in Thoul & Weinberg (1995) and KI. At each time- 
step, the radiative transfer is solved with the method devised 
by TU, in which both absorption and emission (scattering) 
of ionizing photons are explicitly taken into account. We 
assume for simplicity that the baryonic gas is composed of 
pure hydrogen. Neglecting helium causes only a minor ef- 
fect in the ionization degree less than about the order of 
10% (Osterbrock 1989; Nakamoto, Susa & Umemura 1998). 
We further assume ionization equilibrium among photoion- 
ization, collisional ionization, and recombination. This as- 
sumption is well justified in the present analysis (see §2.2 of 
KI for discussion). The number of mass shells is Nh = 200 
for baryonic gas and iVd = 2000 for dark matter. At each 
radial point, angular integration of the radiative transfer 
equation is done over at least 20 bins in 6 — — it, where 6 
is the angle between the light ray and the radial direction. 
This is achieved by handling 300 — 700 impact parameters 
for light rays. The radiation field and the ionization states in 
the cloud interior are solved iteratively until the HI fraction, 
Xhi, in each mesh converges within an accuracy of 1%. 

The external UV field is presumed to be isotropic and 
to have a power-law spectrum: 

J{v) = -hi ( — ) " x 10~ 21 erg s" 1 cm" 2 str" 1 Hz" 1 , (1) 

V ^HI / 

where J21 is the intensity at the Lyman edge of hydrogen 
(hum = 13.6eV) and a is the spectral index. We consider 
two typical cases for a, i.e., a = 1 representative for black 
hole accretion and a — 5 for stellar UV sources. Obser- 
vations of the proximity effect in the Lya forest suggest 
J21 = 10 ±0 ' 5 at z = 1.7 - 4.1 (e.g., Bajtlik, Duncan & Os- 
triker 1988; Bechtold 1994; Giallongo et al. 1996), but its 
value is still uncertain at other redshifts. In what follows we 
give the onset of the UV background to be at zuv = 20 and 
study the following two cases for J21: 

(i) Constant UV 

J21 = 1 z < zuv, (2) 

(ii) Evolving UV 

( exp[— (z — 5)] 5 < z < zuv 
J21 = < 1 3 < z < 5 (3) 




< z < 3. 



The form of the UV evolution at z > 5 in (ii) is roughly con- 
sistent with the results of recent models for the reionization 
of the universe (e.g., Ciardi et al 2000; Umemura, Nakamoto 
& Susa 2000). 

We start the simulations when the overdensity of a 
cloud is still in the linear regime, adopting the initial and 
boundary conditions described in KI. The initial overdensity 
profile is Si(r) — <5i(0) sm(kr)/kr, where k is the comoving 
wave number, and the central overdensity <5i (0) is fixed at 
0.2. The outer boundary is taken at the first minimum of 
5i(r), i.e., kr = 4.4934, within which the volume averaged 
overdensity S(< r) vanishes. Following Haiman, Thoul & 
Loeb (1996), the characteristic mass of a cloud M c i ou d is 



defined as the baryon mass enclosed within the first zero 
of Si(r), i.e., kr — n. Collapse redshift z c is defined as the 
epoch at which M c i ou d would collapse to the center in the 
absence of thermal pressure. Circular velocity V c and virial 
temperature TUr are related to z c and M c i ou d via usual def- 
initions: 

"'= I5 - 9 (^ ! w)"* <1+&) " 1 km '"- (4) 

T*. (jt,) (^g^V.) K,,5) 

where fi is the mean molecular weight in units of the proton 
mass m p . 

The collapse of a gas shell is traced until it reaches the 
rotation radius specified by the dimensionless spin parame- 
ter; 



r I0 t = 




where r ta is the turnaround radius of the gas shell, and we 
adopt a median for the spin parameter, A ta = 0.05 (Efs- 
tathiou & Jones 1979; Barns & Efstathiou 1987; Warren et 
al. 1992). Below the size given by (^|), the system would 
attain rotational balance and forms a disk eventually. 



3 RESULTS 

3.1 Significance of radiative transfer 

To demonstrate the significance of coupling the radiative 
transfer of ionizing photons with hydrodynamics, we first 
present in Fig. 1 estimations for relevant timescales at the 
center of a uniform static cloud with circular velocity V c 
and collapse epoch z c , where the external UV is specified by 
J21 and a spectral index a. The photoionization timescale 
ti on and the photoheating timescale thcat are compared to 
the dynamical timescale tdyn, either by solving the radiative 
transfer of UV photons through the cloud or by just as- 
suming the optically thin medium. Wherever necessary, we 
have adopted the temperature T = 10 4 K, the proton num- 
ber density n H = n£ r = 5.0 x 10- 5 (n b ft 2 /0.025)(l + z c ) 3 
cm -3 , the mass density p = m p riHO,o/0.b, and the radius 
R = R vil = 64(U c /30km s- 1 )(n /i 2 /0.25)- 1/2 (l + z c )" 3/2 
kpc. In the low-density limit, the contours of t[ on — t^yn 
and thcat = idyn both approach asymptotically the optically 
thin case J21 oc -J rig r oc (1 + z c ) 3 ^ 2 . The differences orig- 
inated from the cloud sizes are rather small compared to 
those arisen when one incorporates the radiative transfer or 
not. 

Fig. 1 predicts that the significance of radiative trans- 
fer effects increases with increasing redshift (i.e., increasing 
cloud density) and decreasing J21 for a given a. Under the 
UV evolution given in equation (^), for example, td yn < iion 
and idyn < thcat are satisfied at the center of a cloud with 
V c = 30km s _1 at z c ^ 10, and the cloud is shielded against 
the external radiation in terms of photoionization as well 
as photoheating. At 6 ^ z c 10, thcat < td yn < iion is 
achieved, indicating that the cloud center is heated but not 
ionized by the UV radiation (Gnedin & Ostriker 1997; KI). 
Contrastively, both ii on and iheat become shorter than tdyn 
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at z c < 6 and the cloud is likely to evolve in a similar fashion 
to the optically-thin case. We can see qualitatively the same 
relations for a — 5, except that shielding of ionization and 
heating occurs almost simultaneously at lower redshift. It 
should be noted that the above estimations are made for a 
virialized cloud using the UV intensity only at z c . As shown 
in forthcoming sections, the central density of a collapsing 
cloud actually continues to ascend to above and the ra- 
diative transfer effects can become important even at low 
redshifts. In addition, possible changes of the UV intensity 
during the dynamical growth of a cloud can also affect its 
ionization structure. 



shift, the whole cloud is kept self-shielded in terms of both 
photoionization and photoheating at z ^ 12. For a cloud 
with relatively high z c (z c = 4.8, Fig. ^a), the cloud cen- 
ter begins to contract before the penetration of the external 
UV. Hence, the dynamical evolution closely coincides with 
that without the UV background. As the virial temperature 
of the cloud is T v i r ~ 3 x 10 4 K, the could center is shock- 
heated to above 10 4 K and cools via atomic cooling. On the 
other hand, a cloud with low z c (z c — 0.5, Fig. ^j) is once 
photoionized to the similar level to the constant UV case at 
3 < z < 5, but is able to collapse as the UV intensity drops 
at lower redshifts. 



3.2 Dynamical evolution 

Fig. ^ shows the dynamical evolution of a cloud with 
V c = 29km s _1 under the constant UV background at 
z < zuv = 20. This cloud would collapse at z c = 3 if 
there were no UV background. In practice, the cloud turns 
around and contracts in the inner parts, while it continues 
to expand in the outer envelope. For a hard UV spectrum 
with (J2i,a) = (1, 1), the gas is ionized and heated up to 
T ~ 10 4 K promptly at the onset of the UV background. 
For a soft UV spectrum with (Jai,a) = (1,5), the cloud 
center is kept self-shielded against the external field and the 
temperature ascends more gradually. 

For comparison, results of the optically-thin and pure 
absorption calculations are also presented in Fig. ^. The 
pure absorption case is based on the analytical formalism 
described in KI. The cloud evolution is altered in no small 
way by different treatments of the UV radiation. Under the 
optically-thin assumption, the cloud is completely prohib- 
ited from collapsing for (J2i,a) = (1,1), and from being 
self-shielded for (J2i,a) = (1,5). The pure absorption ap- 
proximation leads to accelerating the collapse of the central 
core and to underestimating photoionization and photoheat- 

illg ' 

Fig. H exhibits the radial profiles of the same clouds at 
z = 3 for (Jai,a) = (1,5). The inner part of a cloud does 
not settle into hydrostatic equilibrium, but rather undergoes 
isothermal run-away collapse and the density profile follows 
the self-similar solution of Bertschinger (1985). The trend is 
regardless of the treatment of the UV radiation transfer, be- 
cause the temperature of thermal equilibrium in optically- 
thin media comes close to 10 4 K in the high density limit. 
The present results demonstrate that the approximation of 
hydrostatic equilibrium employed in previous analyses (Kep- 
ner, Babul & Spergel 1997; Barkana & Loeb 1999) breaks 
down for the final states. 

The ionization structure is quite different depending on 
whether one incorporates the radiative transfer or not. In an 
optically-thin cloud, the ionization degree changes merely 
according to the ionization parameter, ns/J2i, whereas the 
radiative transfer effects produce a self-shielded neutral core 
(Fig. |^c). Fig. ^(d) further indicates that the UV heating 
rate is reduced significantly in the self-shielded region, with 
an increasing contribution of scattered photons to the total 
heating rate. Implications of the present results on H2 cool- 
ing [thin lines in panel (d)l will be discussed in detail in Sec 

Effects of the evolution of the UV background are illus- 
trated in Fig. M. Due to very low UV intensity at high red- 



3.3 Criteria for collapse 

Figs U and ^ summarize the results of present calculations 
for a variety of initial conditions on a 14 — z c plane. As our 
simulations assume ionization equilibrium and only incor- 
porate atomic cooling, they can be most securely applied to 
a cloud once heated to above 10 4 K either by shock or by 
UV photons during the course of its evolution. In contrast, 
evolution of lower temperature systems is still tentative and 
may be altered once cooling by molecular hydrogen is explic- 
itly taken into account. These figures thus distinguish "high 
temperature clouds" (circles) defined as those photoheated 
to > 10 4 K or those with T vil > 10 4 K, and the other "low 
temperature clouds" (triangles). Each of these populations 
are further classified by open and filled symbols depend- 
ing on whether or not they collapse to the rotation barrier 
within the present age of the universe. 

It is obvious that the UV background prohibits small 
clouds from collapsing even if the transfer of ionizing pho- 
tons is considered. Under a constant UV flux (Fig. ^), 
threshold circular velocity for the collapse gradually in- 
creases with decreasing redshift, i.e., decreasing cloud den- 
sity. The threshold velocity also decreases with increasing 
photon spectral index a because of the smaller number of 
high energy photons. At z c ^ 5, the threshold falls below 
10 4 K, because the cloud center starts to contract before 
the onset of the UV background and remains impervious to 
the external photons. In the presence of the UV evolution 
(Fig. ^), the threshold velocity increases sharply at z c £ 3 
and drops slightly at z c ,$ 1. The central temperature of a 
cloud denoted by an open triangle can reach ~ 10 4 K by adi- 
abatic compression and the collapse is promoted by atomic 
cooling. A cloud denoted by a filled triangle fails to collapse 
because of the lack of coolant at T < 10 4 K in our simu- 
lations. As will be discussed in Sec ^, however, evolution 
of these "low temperature clouds" may be modified by H2 
cooling. 

KI have suggested that the threshold for collapse is 
roughly determined by the balance between the gravita- 
tional force and the thermal pressure gradient when the gas 
is maximally exposed to the external UV flux. To confirm 
this, we plot in the same figures a relation T v i r = T,iq ax , 
where Tcq ax is the maximum equality temperature defined 
semi-analytically as follows. Firstly, given the initial over- 
density profile, the collapse of a spherical perturbation is 
approximated by the self-similar solution of Bertschinger 
(1985). Secondly, at an arbitrary stage of the collapse, one 
can compute the equality temperature at which radiative 
cooling balances photoheating in the optically thin limit, 
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using the central density deduced from the self-similar so- 
lution and the UV intensity assumed in the simulation. Fi- 
nally, T<^ ax is set equal to the maximum of such tempera- 
ture. For the gas exposed to a constant UV flux from the 
linear regime, T™ ax is attained essentially at turn-around. If 
a cloud has T v i r above T™ ax , it is likely to be gravitationally 
unstable against gas pressure. 

Figs H and M show that the above criteria agree rea- 
sonably well with our simulation results based on the ac- 
curate treatment of the radiative transfer. This is because 
the dynamical evolution is basically regulated by the Jeans 
criterion when a cloud is heated up to ~ 10 4 K by the UV 
background before contraction. At high redshifts (z c £ 6 in 
Fig. H and z c Si 3 in Fig. |^), however, the cloud center is 
strongly self-shielded from early stages and the approxima- 
tion of the optically thin limit breaks down. The Jeans scale 
at these epochs is reduced below the relation T v i r = T™ ax . 



4 IMPLICATIONS FOR GALAXY 

FORMATION IN THE UV BACKGROUND 

The present RHD calculations give an accurate prediction 
for the suppression of pregalactic collapse due to the UV 
background, which has been one of the primary concerns 
from the viewpoint of galaxy formation (Umemura & Ikeuchi 
1984; Ikeuchi 1986; Rees 1986; Bond, Szalay & Silk 1988; Ef- 
stathiou 1992; Babul & Rees 1992; Zhang, Anninos & Nor- 
man 1995; Thoul & Weinberg 1996; KI). What is of addi- 
tional significance in this context is the subsequent forma- 
tion of stars in collapsing clouds. In order for stars to form, 
a cloud needs to be cooled down to below 10 4 K by hydro- 
gen molecules, because they are the only coolant in metal- 
deficient gas (e.g., Peebles & Dicke 1968; Matsuda, Sato & 
Takeda 1969; Tegmark et al. 1997). H2 cooling is a two- 
body collision process (e.g., Hollenbach & Mckee 1989; Galli 
& Palla 1998), while the photoheating rate is in proportion 
to the density. The potentiality of H2 cooling is therefore 
enhanced with increasing density. In this respect, run-away 
collapse should provide favorable situations for H2 cooling. 

Based on our simulation results presented in the pre- 
vious section, we further investigate the possibility of star 
formation in a collapsing core. We suppose that cooling be- 
low 10 4 K becomes efficient if the following two conditions 
are both satisfied; 1) photodissociation of molecular hydro- 
gen by the UV photons in the Lyman- Werner bands at 
11.26-13.6keV (e.g., Stecher & Williams 1967) is less impor- 
tant than other H2 destruction processes, and 2) H2 cooling 
overtakes UV photoheating. The first condition is fulfilled 
for clouds heated up to > 10 4 K, because destruction of H2 is 
dominated by collisional dissociation insofar as T <; 2000K 
(Corbelli, Galli & Palla 1997). More specifically, requiring 
that the timescale of H2 dissociation via collisions with H + 
(reaction 12 of Shapiro & Kang 1987) is shorter than that of 
photodissociation (Draine & Bertoldi 1996; Omukai & Nishi 
1999) in the conservative optically thin limit, we have for 
the electron density 

n c > 4.7 x 1(T 5 f ^ff ) J21 exp(21200 K/T) cm -3 , (7) 

which is amply satisfied at the collapsing core of our simu- 
lated clouds. The second condition depends on the amount 



of H2 formed. In the absence of the external UV field, 
the H2 abundance in the metal-deficient postshock layer 
converges roughly to Xh2 ~ 1CP 3 (e.g., Shapiro & Kang 
1987; Ferrara 1998; Susa et al. 1998). Under the UV field, 
Xn2 ~ 1CP 3 is also achieved if photoheating is strongly at- 
tenuated by self-shielding, while the abundance is reduced 
down to Xh2 ~ 10 -6 in the case of weak attenuation of 
photoheating (Kang & Shapiro 1992; Corbelli et al. 1997; 
Susa & Umemura 2000; Susa & Kitayama 2000). It is also 
likely that the H2 abundance depends on the spectrum of 
the impinging radiation and the details of the ionizing struc- 
ture of the medium (e.g., Kang & Shapiro 1992; Ciardi et al. 
2000). Leaving such complexities elsewhere (Kitayama et al. 
in preparation), we make crude estimations for the cooling 
rate using the H2 cooling function of Hollenbach & Mckee 
(1989) and assuming A"h2 = 10 -3 in the full transfer and 
pure absorption cases, and Xn2 = 10 -6 in the optically thin 
case. 

Fig. ^(d) shows that H2 cooling greatly overwhelms pho- 
toheating in the self-shielded neutral core for Xu2 = 10 -3 . 
It turns out that what enables H2 cooling to be effective is 
essentially the attenuation of photoheating rather than the 
final abundance of H2 molecules. In fact, in the strongly self- 
shielded core, H2 cooling remains to be dominant even if the 
abundance is reduced to a level similar to the optically-thin 
case, i.e., Ah2 = 10 -6 (the H2 cooling rate is roughly propor- 
tional to Xh2 at T ^ 10 4 K). We have further checked that 
this is the case with every collapsed cloud plotted as an open 
circle in Figs ^ and [| Under the optically-thin approxima- 
tion, on the other hand, UV heating is much stronger than 
H2 cooling, prohibiting the cooling below 10 4 K. 

While the above arguments support the possibility of 
H2 cooling during the collapse of "high temperature clouds" 
(open circles in Figs ^ and g) , the situation is rather intricate 
for "low temperature clouds" (triangles in the same figures) . 
Once incorporated formation and destruction of molecular 
hydrogen explicitly in our calculations, these clouds would 
be able to collapse as long as the external UV flux is negligi- 
ble or very weak. For somewhat stronger UV flux, photodis- 
sociation may operate efficiently to suppress H2 formation 
and to prohibit the collapse (Haiman, Rees & Loeb 1996, 
1997; Haiman, Abel & Rees 2000; Ciardi et al. 2000). Al- 
ternatively, self-shielding induced by run-away collapse may 
still enable a sufficient amount of H2 to be formed for the 
system to cool down to ~ 100K. These points will be inves- 
tigated thoroughly in our future work (Kitayama et al. in 
preparation) . 

In summary, self-shielding against the UV background 
is indispensable for H2 cooling and subsequent star for- 
mation to proceed. Run-away collapse is likely to promote 
clouds with T v i r > 10 4 K to cool down before collapsing to 
the rotation barrier given by equation (^|). These results are 
closely relevant to the formation of dwarf galaxies at high 
redshifts as well as to metal injection into intergalactic space 
(e.g., Nishi & Susa 1999). Note that the above conclusions 
do not conflict with those of several recent works (Haiman 
Rees & Loeb 1996, 1997; Ciardi et al. 2000; Haiman et al. 
2000) that the so called "radiative feedback" is operating 
in the early universe and suppresses the collapse of some 
objects via photodissociation of H2 molecules. Current sim- 
ulations are mainly aimed to study the clouds once heated 
to above 10 4 K either by shock or by UV photons, for which 
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photodissociation has minor impacts on H2 formation com- 
pared to collisional dissociation (eq. 0). In addition, run- 
away collapse yields a highly self-shielded core which is par- 
ticularly favorable for H2 formation. Photodissociation may 
be of greater significance in smaller clouds and lower density 
regions. 

Finally, we discuss the effects of the geometry of cloud 
evolution upon H2 cooling. Near and above the Jeans scale, 
clouds are expected to contract spherically to a first order 
approximation and result in run-away collapse as shown in 
the present paper. Clouds far above the Jeans scale, on the 
other hand, are likely to undergo pancake collapse. The pan- 
cakes would not end up with run-away collapse, because the 
gravity of the sheet is readily overwhelmed by the pressure 
force. Susa & Umemura (2000) have explored pregalactic 
pancake collapse including the UV transfer and H2 forma- 
tion, and shown that the collapsing pancakes bifurcate, de- 
pending upon the initial masses, into less massive UV-heated 
ones and more massive cooled ones. This is because the self- 
shielding in sheet collapse is governed by the column density 
of a pregalactic cloud. To conclude, star formation in pre- 
galactic clouds in the UV background is strongly regulated 
by the behaviour of collapse and the manner of the radiation 
transfer. 
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Figure 1. J21 — z c diagram on timescales relevant to the evolution of a prcgalactic cloud in the UV background with (a) a = 1 and 
(b) a = 5. The photoionization time t[ on , the photohcating time theati an< i the dynamical time t^ yn are evaluated at the center of a 
stationary uniform cloud with V c and z c as described in the text. Lines indicate the contours tj on = td yn (dotted) and ii leat = tdyn 
(solid) for V c = 100 km s — 1 (thick), 30 km s — 1 (mildly thick), and in the optically thin limit (thin). In the region under these contours, 
t dyn becomes shorter than t; on and ti loat , respectively. For reference, the track of the UV evolution given in equation (jj) is marked by 
open circles. 
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Figure 2. Evolution of radius (top panels), temperature (middle panels), and HI fraction (bottom panels) of gas shells enclosing 0.5% 
(thick lines) and 90% (thin lines) of total gas mass in a cloud with V c = 29 km/s and z c = 3 under the constant UV background; (a) 
(J21, a) = (1, 1) and (b) (1, 5). Different line types indicate the full transfer (solid), optically-thin (dotted), and pure absorption (dashed) 
cases, respectively. Also plotted in the top panels are the no UV (dot-dashed) results. 
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Figure 3. Radial profiles at z = 3 of (a) hydrogen density, (b) temperature, (c) HI fraction, and (d) UV heating rate, for (J21, a) = (1, 5), 
V c = 29 km/s, and z c = 3. Different line types indicate the full transfer (solid), optically- thin (dotted), and pure absorption (dashed) 
cases, respectively. The heating rate caused by scattering in the full transfer case is added to panel (d). Also plotted in (d) are estimations 
for the H2 cooling rate inside the shock front at r ~ lOkpc, assuming A"h2 = 10 — 3 in the full transfer or pure absorption cases, and 
Xjj2 = 10~ 6 in the optically-thin case. 
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Figure 4. Similar to Fig. |2| except that the effects of the UV evolution (solid lines) are illustrated in comparison to the constant UV 
(dotted) and no UV (dot-dashed) cases; (a) z c = 4.8, V c = 29 km/s, and (b) z c = 0.5, V c = 29 km/s. Full transfer results are plotted 
assuming a = 1. 
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Figure 5. V c — z c diagram describing collapse and self-shielding from our full-transfer RHD calculations; (a) (J2i,a) = (1, 1), and (b) 
(1,5). Circles are the present secure results for "high temperature clouds" and triangles are the tentative results for "low temperature 
clouds" (see text for definitions). Filled circles are clouds which are prohibited to collapse due to the UV heating, and open circles are 
those which undergo run-away collapse with T v ; r > 10 4 K. Filled triangles are clouds which cannot collapse to the rotation barrier, and 
open triangles are those which collapse due to atomic cooling. Also plotted by dashed lines are the relation T v ; r = T™ ax defined in the 
text. 
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Figure 6. Same as Fig. |s| except that the UV evolution is taken into account. 
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